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ON  THE  BUCKLING  FORCE 
OF  FLOATING  ICE  PLATES 


Arnold  D.  Kerr 


INTRODUCTION 

The  analytical  determination  of  the  largest  horizon¬ 
tal  force  a  floating  ice  cover  may  exert  on  a  structure 
is  based  on  the  criterion  that  an  ice  force  cannot  be 
larger  than  the  force  capable  of  breaking  up  the  ice 
cover. 

In  many  publications,  as  discussed  by  Korzhavin 
(1971)  and  Michel  (1970),  this  failure  criterion  is  re¬ 
lated  to  a  “crushing  strength”  of  the  ice  cover.  There 
it  is  assumed  that  the  failure  mechanism  consists  of  the 
crushing  and/or  splitting  of  the  floating  ice  plate  in  the 
immediate  vicinity  of  the  structure. 

However,  in  a  number  of  laboratory  and  field  tests 
it  was  observed  that  for  ice  covers  that  are  relatively 
thin  (compared  to  the  width  of  contact  of  plate  and 
structure)  the  ice  cover  failed  by  buckling  in  the  vicin¬ 
ity  of  the  contact  area.  Thus,  for  these  cases,  the  buck¬ 
ling  force  is  smaller  than  the  crushing  force,  and  the 
determination  of  the  largest  force  should  be  based  on 
the  force  at  which  the  floating  ice  cover  is  expected  to 
buckle. 

The  related  buckling  analyses  and  tests  were  recently 
reviewed  and  discussed  by  Kerr  (1978).  In  all  these 
analyses  it  was  assumed  that  the  buckling  load  corres¬ 
ponds  to  the  first  bifurcation  load  pct,  an  assumption 
justified  when  the  post-buckling  equilibrium  branch 
rises  monotonically  with  increasing  displacements. 
However,  according  to  recent  studies  by  Lekkerkerker 
(1962),  Kerr  (1972),  and  Plaut  (1978),  this  is  not  the 
case  in  general.  The  purpose  of  the  present  report  is  to 
clarify  this  situation  for  the  floating  ice  plate  problem. 


SIMPLE  MODEL  FOR  ICE  COVER  BUCKLING 

To  study  the  basic  features  of  die  title  problem  we 
consider  first  the  simple  model  shown  in  Figure  1,  which 


consists  of  “rigid  bars”  interconnected  by  elastic  spiral 
springs  and  resting  on  straight  vertical  springs.  This 
model  contains  the  observed  buckling  mechanism  of  a 
compressed  floating  ice  cover  and  is  amenable  to  a 
simple  nonlinear  exact  analysis.  It  is  assumed  that  the 
rigid  bars  are  horizontal  when  subjected  to  the  uniform 
load  distribution  q0  (their  own  weight).  In  this  position 
of  equilibrium  the  spiral  springs  are  stress-free  and  each 
of  the  straight  vertical  springs  is  compressed  by  a  force 
kw0.  Thus,  the  vertical  displacement  is 


w0=q0Uk  (1) 

where  k  is  the  stiffness  of  the  vertical  springs.  Next, 
the  bars  are  subjected  to  an  axial  force  P.  For  increas¬ 
ing  P,  deformed  states  of  equilibrium  may  exist,  as 
shown  in  Figure  2a.  Assuming  that  the  bars  at  joint  2 
will  not  separate  from  the  base  spring  and  thus  that 

w0jJZ.sinO  or  qo/k£sin0  (2) 

the  equilibrium  equation  for  the  deformed  state  is 

P*  sin0  =  d+k*sir\0  cos0  (3) 

where,  denoting  the  parameter  of  the  spiral  springs  by  s, 

P*  =  PL/3s ;  k*  =  kl?l6s.  (4) 

Note  that  q0  does  not  appear  in  equilibrium  equation  3 
because  of  the  linearity  of  the  base  springs.  Also  note 
that  eq  3  is  satisfied  for  0  =  0.  Thus,  the  straight  state 
is  in  equilibrium  (but  not  necessarily  stable)  for  any 
load  P.  For  very  small  values  of  0,  eq  3  reduces  to 

(/>*_1_A*)0  =  O.  (3') 

Thus,  at  the  bifurcation  point 


L 


\ 


j 


undeformed  state (q0*0) 


Figure  1.  Simple  buckling  model  (q0  is  the  uniform  load  distribution, 
s  the  spiral  spring  parameter,  and  L  the  distance  between  the  vertical 
springs). 


(b)with  lift- off 


P*x^Uk*.  (5) 

Equation  3  was  evaluated  numerically  for  various  val¬ 
ues  of  k*.  The  results  are  shown  in  Figure  3. 

For  a  proper  interpretation  of  the  results  shown  in 
Figure  3a,  note  that,  as  shown  in  Figure  4  (Kerr  1970), 
the  equilibrium  branch  rises  monotonically  when  the 
model  bars  are  constrained  only  by  a  spiral  spring 
(bending  effect),  whereas  the  equilibrium  branch 
drops  when  the  model  bars  are  constrained  by  a 
straight  spring  (effect  of  base).  This  is  the  reason  why 
in  Figure  3a,  for  large  values  of  k*  (i.e.  when  the  base 
stiffness  predominates),  the  equilibrium  branch  drops 
at  first  and  then  rises,  exhibiting  a  lower  buckling  load 

Pl<Pc  r 

From  the  graphs  in  Figure  3a,  it  follows  that  the 
value  of  k*  which  separates  the  different  responses  is 
located  in  the  interval  0  <  k*  <  0.5.  Using  the  per¬ 
turbation  method,  Kerr  (1972)  showed  that  for  the 
model  under  consideration 

*?*'/*.  (6) 


Consider  the  model  with  k  *  =  1 .0.  According  to 
Figure  3,  for  each  P  <  PL  there  exists  only  one  straight 
equilibrium  state,  for  PL<  P<  Pcr  there  exist  five 
equilibrium  states,  and  for  each  P  >  Per  three  states  of 
equilibrium  are  possible.  It  may  be  shown  (Kerr  1974) 
that  a  straight  equilibrium  state  is  stable  below  PCT  and 
unstable  above  it,  that  the  equilibrium  states  on  branch 
PcrL  are  unstable,  and  that  the  equilibrium  states  on 
branch  LB  are  stable. 

When  the  axial  load  P  has  a  small  vertical  eccentric¬ 
ity,  the  equilibrium  branches  are  as  shown  by  the 
dashed  lines  in  Figure  3b.  The  effect  of  vertical  shape 
imperfections  is  similar.  Note  that  with  increasing 
small  eccentricities,  or  shape  imperfections,  the  upper 
buckling  load  Pu  drops  rapidly.  Structures  of  this  type 
are  referred  to  in  the  literature  as  "imperfection  sensi¬ 
tive.” 

Because  of  the  existence  of  load  and  geometrical 
imperfections  (and  dynamic  inputs)  in  an  actual  situa¬ 
tion,  buckling  will  take  place  for  PL  <  P  <  PCT ,  where 
/>[_  is  the  value  for  the  structure  with  imperfections. 
Thus,  for  the  cases  that  exhibit  a  PL  value,  Pcr  is  not 
the  true  buckling  load. 
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Figure  3.  Equilibrium  branches  for  model  without  lift-off.  For  0  <  0  the  branches  are 
symmetrical  to  the  P*  axis. 


P  P 


Figure  4.  Simple  models  and  their  equilibrium  branches. 


Next  consider  the  case  when  the  bars  are  resting  on,  The  reason  for  the  large  post-buckling  deformations  in 

but  are  not  attached  to,  the  base  spring  at  joint  2;  i.e.  Figures  3  and  5  is  the  assumption  that  an  axial  force  P 

they  may  lift  off  the  spring  base,  as  shown  in  Figure  remains  the  same  also  at  the  post-buckling  equilibrium 

2b.  This  may  occur  in  an  actual  situation  with  floating  points.  Flowever,  it  was  shown  by  Kerr  (T973)  that, 
ice  covers.  For  this  case,  until  lift-off  takes  place,  the  when  the  axial  compression  force  is  induced  by  a  tern- 

governing  equilibrium  equation  is  eq  3.  After  lift-off,  perature  raise,  the  axial  force  drops  and  the  post-buck  - 

namely  when  ling  deformations  are  much  smaller. 

According  to  Figure  5a  the  admissibility  of  lift-off 
w0  <  L  sin0  or  q0/k  <  sinfi  (7)  substantially  reduces  the  PL  value,  especially  for  small 

values  of  q*.  The  resulting  equilibrium  branch  for  q* 

the  equilibrium  equation  is  =  0.1  is  shown  as  a  solid  line.  Note  that  the  correspond¬ 

ing  equals  1 .2;  thus  this  value  is  about  25%  lower  than 
P*  sin#  =  6  +  q*  cosO  (8)  the  P£  for  the  case  when  lift-off  is  not  allowed. 

Figure  5b  shows  that  the  admissibility  of  lift-off  cre- 

where  ates  a  PL  smaller  than  Pcr ,  even  for  an  equilibrium 

branch  that  is  monotonically  increasing  with  no  lift-off. 
P*  =  PL/3s;  q*  =  q0L2/6s.  (9)  This  is  an  example  where  a  stability  analysis  of  the  bi¬ 

furcation  point  on  the  undeformed  branch  (for  example, 
The  corresponding  equilibrium  branches  for  k  *  =  0-2,  the  Koiter  method)  is  not  suitable  for  predicting  whether 

1 .0  and  a  range  of  q*  values  are  shown  in  Figure  5.  PL  exists. 
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Figure  5.  Equilibrium  branches  for  mode!  with  lift-off . 


The  findings  of  this  section  suggest  that  for  com¬ 
pressed  floating  ice  plates,  even  in  the  absence  of  lift¬ 
off  between  plate  and  liquid  base,  apL  smaller  than 
pcr  may  exist,  as  shown  in  Figure  3a.  The  results  in 
Figure  5  suggest  that  the  possibility  of  lift-off  may 
greatly  affect  the  post-buckling  equilibrium  branches 
and  thus  also  the  lower  buckling  loads pL.  Occurrence 
of  lift-off  may  even  create  a  pL  value  where  one  does 
not  exist  when  lift-off  is  not  allowed. 

To  establish  whether  floating  ice  plates  do  respond 
in  this  way  requires  the  solution  of  the  nonlinear  equa¬ 
tions  for  plates  on  a  liquid  base.  The  study  of  all  these 
phenomena  is  beyond  the  scope  of  the  present  report. 
Therefore,  in  the  remainder  of  this  report  only  the  be¬ 
havior  of  the  equilibrium  branches  of  compressed 
floating  ice  plates  in  the  vicinity  of  pcr  (as  shown  in 
Figure  3a)  is  investigated  in  order  to  establish  if,  and 
under  what  conditions,  a pL  smaller  thanpcr  will  exist 
in  the  absence  of  lift-off. 


SEMI-INFINITE  FLOATING  PLATE 
WITH  FREE  EDGE 

The  problem  under  consideration  is  shown  in  Figure 
6.  Because  the  load  p  (per  unit  length  of  plate  boundary) 
is  assumed  to  be  constant,  the  plate  displacements  will 
be  functions  of  x  only.  The  vertical  displacement  com¬ 
ponent  of  a  point  (x,  y)  on  the  middle  plane  is  assumed 
in  the  form  w0  +  w(x);  w0  =  constant  is  the  rigid  body 
displacement  due  to  the  weight  of  plate  q0  (per  unit 
area),  and  w(x)  is  caused  by  bending  deformations. 

Thus,  q0  =  yw0.  Following  the  derivations  in  Kappus 
(1939)  or  Kerr  (1972),  using  Lagrange  coordinates  and 
neglecting  the  extensibility  of  the  middle  plane,  the 
differential  equation  for  w(x)  is 

ntt  a  *  n  ,fr  i  "3  a  '2  "3 

w  ww  w  +w  4.4 w  w 

(1  -w'2)2  (1-w'J)* 


+  £ 


D  (1  -w'*  )3,i 


+  2  w  =  0 
D 


(10) 


Figure  6.  Semi-infinite  plate  with  free  edge. 


where  D  =  Eh3 1 12(1  -u2)  is  the  flexural  rigidity  of  the 
plate,  7  is  the  base  modulus,  and  ()'  =  d(  )/dx.  The 
boundary  conditions  are 


w"(0)  =  0 

(ii) 

(12) 

The  formulation  is  highly  nonlinear.  In  the  mathe¬ 
matical  literature  there  are  no  methods,  as  yet,  for 
obtaining  an  exact  closed  form  solution  for  this  bound¬ 
ary  value  problem.  Therefore,  in  the  following,  we 
analyze  the  post-buckling  equilibrium  branch  in  the 
vicinity  ofpcr  in  order  to  establish  if  pL  exists.  This 
is  done  using  the  perturbation  method  (Kerr  1972). 

Following  Keller  (1968)  the  variablesp  and  w(x) 
are  expanded  in  a  Taylor  series: 

p(c)=p(0)+pe(0)e+pee(0)  e2/2!  +pete( 0)  e3/3!  +... 

w (x,  e)  =  w (x,  0)  +  we (x,  0)e  +  w(( (x,  0)  e1  /2! 

+  vveff(x,0)  e3/3!  +...  (13) 

where  e  is  a  small  parameter.  Denoting  symbolically 
the  nonlinear  formulation  consisting  of  eq  1 0-1 2  by 

F[w(x),p}=  0,  (14) 

it  then  follows,  because  of  eq  1 3,  that  F  =  F(e)  only. 
Assuming  that  F  is  differentiable,  we  write 

F(e)  =  F (0)  +  F((0)e  +  Fee{0)  e2/2!  +...  =  0. 

05) 

Since  e  is  an  arbitrary  parameter,  the  above  equation 
is  satisfied  when 

F(0)  =  0;  Ft  (0)  =  0;  Fj  0)=0;...  (16) 

These  are  the  equations  for  the  determination  of  the 
coefficients  in  eq  1 3.* 

It  should  be  noted  that  in  eq  1 3  the  term  w(x,  0) 
represents  the  pre-buckling  state.  Hence 

w{x,  0)=0.  (17) 

The  formulation  which  corresponds  to  the  first  con¬ 
dition  in  eq  16,  F( 0)  =  0,  is  satisfied  for  any  point  on 
the  equilibrium  branch  with  w(x,  0)  =  0. 


Applying  the  second  condition  in  eq  16,  Fe( 0)  =  0, 
to  the  eq  10-1 2  and  noting  that  because  of  eq  1 7  the 
derivatives  (x,  0)  =  0,  we  obtain 


Lwe(x,  0)=0  0^x<°° 

(18) 

where 

Z.=D()iv+p(0)()"+7 

(19) 

and  the  boundary  conditions 

o 

It 

o' 

cT 

V 

(20) 

Dw'J'  (0, 0)  +p(0)w£  (0, 0)  =  0. 

(21) 

This  is  a  linear  eigenvalue  problem  for  we(x,  0).  It  is 
identical  to  the  formulation  in  the  Euler  method  for  ob 
tainingpcr. 

This  eigenvalue  problem  is  satisfied  by  the  trivial 
solution 

we(x,0)  =0 

(22) 

which  represents  the  continuation  of  the  undeformed 
equilibrium  branch.  It  is  also  satisfied,  for  the  equili¬ 
brium  branch  which  corresponds  to  w(x,  e)  +  0,  by 
the  first  eigenvalue  (Rzhanitsyn  1955): 

p{0)=pcr=yJyD 

(23) 

and  the  corresponding  eigenfunction: 

we(x,  0)  =  fle^[sin(\/3ax)--\/3cos(\/3ax)] 

(24) 

where  B  is  an  arbitrary  constant  and  a  =  tyy)l6b. 

Next,  we  form  Fee  ( 0)  =  0.  The  resulting  equations 
are  the  differential  equation  for  we£(x,0) 

Lwee(x,  0)  =  -2pe(0)w/'(x,0)  0  ^  x  <  °°  (25) 

and  the  boundary  conditions 

<(0,0)  =  0  (26) 

Dw'”  (0, 0)  +p(0)< (0, 0)  =  -2 pe(0)w'((0, 0) 

(27) 

where  the  differential  operator  L  is  given  in  eq  19.  The 
above  boundary  value  problem  is  satisfied  for 


’Equation  15  may  also  be  obtained  by  substituting  eq  13  into 
14  and  then  by  grouping  terms  of  equal  powers  in  t. 


we(x,0)  =  0  and  wee(x,  0)  =  0 


(28) 


which  represent  the  continuation  of  the  undeformed 
equilibrium  branch.  For  the  equilibrium  branch  at 
the  bifurcation  point  which  corresponds  to  w(x,  e)  0, 
the  operator  L  in  eq  19  is  singular  [note  that  the  cor¬ 
responding  homogeneous  problem  is  identical  to  the 
one  for  we(x, 0)].  Hence,  a  solution  we<f(x,0)  will 
exist  only  if  the  Fredholm  alternative  is  satisfied. 

In  this  connection,  the  nonhomogeneity  in  the  se¬ 
cond  boundary  condition  is  eliminated  first,  by  intro¬ 
ducing  a  new  variable  ^  (x)  as  follows: 

%(0,0).  (29) 

P(O) 

With  this  transformation,  differential  equation  25 
becomes 

L  v,  =  2pe(0)jJ^}  *£(0,0) -w'e'(x, 0  J  (25') 

and  the  boundary  conditions  in  eq  26  and  27  become 


vj'(0)  =  0 

(26') 

Dv;"(o)+p(o)v;(o)  =  0. 

(27') 

According  to  the  Fredholm  alternative,  a  solution 
V]  (x),  and  hence  we<,(x,0),  will  exist  only  if 


(30) 


where  \p(x)  is  the  nonzero  solution  of  the  homogeneous 
adjoint  problem. 

Because  the  operator  is  self-adjoint,  it  follows  that 
*P(x)  =  we(x,0)  (31) 

and  the  solvability  condition  (eq  30)  reduces  to 


Pe(0)  we(°.°)  -w''(x,0)]  we(x,0)  dx  =  0. 

(32) 

Noting  that  we(x,0)  is  given  in  eq  24  and  performing 
the  integrations,  it  may  be  shown  that  the  integral 
does  not  vanish.  Hence  eq  32  is  satisfied  when 

Pe(0)=0.  (33) 


given  by  eq  25-27  reduces  to  the  problem  for  wf(x,0) 
given  by  eq  18-21 .  Thus  wee(x,  0)  and  we(x,0)  differ 
by  an  as  yet  undetermined  coefficient  /J,  namely, 

wfe(x,0)  =  0we(x,O).  (34) 

To  determine  the  next  coefficients  in  eq  1 3  we  form 
Feef  (0)  =  0.  The  resulting  equations  are  the  differential 
equation 

Lweee(x,0)  =f(x)  0^x<~  (35) 


where,  noting  thatpe(0)  =  0, 

f(x)  =  -[60(4w'w”w"'  +  w''3  +wf'2wiv) 


+  9pw/2w"  +3peew"]e=0 

(36) 

and  the  boundary  conditions 

%«(o,o)  =  o 

(37) 

OC  (0,0)  +p(0)w;ff(0,0)  =  -3jpee(0)w;(0,0) 

+  D  [  w'2  (0, 0)w"'  (0, 0)  +  2  w'e  (0, 0)  w"2  (0, 0)  ]  j . 

(38) 

The  above  equations  are  satisfied  for 

ws(x)  =0;  wee(x,0)=0;  weee(x,0)  =  0  (39) 

which  represent  the  continuation  of  the  undeformed 
equilibrium  branch. 

For  the  deformed  branch  at  the  bifurcation  point, 
the  above  equations  constitute  a  nonhomogeneous 
boundary  value  problem  with  a  singular  operator.  To 
facilitate  the  use  of  the  Fredholm  alternative,  the  non¬ 
homogeneity  in  the  second  boundary  condition  is 
eliminated  by  introducing  a  new  variable  v(x)  as 
follows: 

w„f(x,0)=v(x)-%^  (40) 

p(0) 

where  {. . .}  is  from  eq  38.  With  this  transformation, 
differential  equation  35  becomes 

Lv  =  f(x)  +  &  [pee(0K(0,0) 


With  eq  33,  the  boundary  value  problem  for  wet  (x,  0) 
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+  Dw?  (0, 0)<'(0, 0)  +  2  Dw’(  (0,  OX2  (0, 0)  ] 

(35') 


and  the  boundary  conditions  of  eq  37  and  38  become 


p"(0)=0 

o  r> 

Dv"'(0)  +p(Q)v'(0)  =  0. 

(39') 

According  to  the  Fredholm  alternative,  a  solution  v{x), 
and  hence  wfef (x,  0),  will  exist  only  when 

/  |f(x)  +  5gj|pff(0)w;(0,0) +Dwf  (0,0)w;"(0,0) 

+  20w;(0,0)w;j(0.0)j|«(x)  dx  —  0  (41) 

where 

<p(x)  =  w((x,0)  (31') 

is  the  nonzero  solution  of  the  homogeneous  adjoint 
problem. 

Performing  the  integrations  indicated  in  eq  41 , 
noting  that  f{x)  is  given  in  eq  36  and  wf(x,  0)  in  eq  24, 
we  obtain 

=  (42) 

where  B  is  an  as  yet  undetermined  coefficient.  Thus, 
according  to  eq  13, 

P(e)  =pcr--^-ifi2e2.  (43) 


CONCLUSIONS 

Equation  43  proves  analytically,  as  predicted  in 
Simple  model  for  ice  cover  buckling,  that  the  post- 
buckling  branch  of  the  semi-infinite  plate  with  a  “free" 
compressed  edge  exhibits  values  of  p  smaller  than  those 
of  pa,  even  in  the  absence  of  lift-off.  According  to  re¬ 
cent  results  by  Plaut  (1978)  this  may  also  occur  for  a 
floating  plate  strip  "simply  supported”  and  compressed 
along  the  straight  parallel  boundaries. 

These  findings  indicate  that  for  floating  ice  covers 
the  critical  buckling  value pcr  should  be  used  with 
caution.  In  this  connection  note  that  the  analysis  pre¬ 
sented  in  Semi-Infinite  flouting  plate  with  free  edge  de¬ 
termines  only  the  behavior  of  the  post-buckling  equi¬ 
librium  branch  in  the  vicinity  of  pa.  The  determina¬ 
tion  of  the  value  pL,  when  one  exists,  requires  the 
solution  of  the  nonlinear  formulation  based  on  eq  10. 

During  buckling  of  ice  covers  lift-off  may  take  place. 
Therefore,  the  analysis  for  the  determination  of  pL 


should  take  this  possibility  into  consideration.  As 
shown  in  the  model  study,  lift-off  may  strongly  affect 
thepL  value  and  hence  also  the  buckling  load. 

The  parameter  k*  in  the  model  shown  in  Figure  2 
represents  the  ratio  of  the  stiffness  of  the  base  to  the 
stiffness  of  the  floating  ice  cover.  For  floating  covers 
the  stiffness  of  the  liquid  base  is  constant,  whereas 
the  plate  stiffness  reduces  with  decreasing  thickness 
of  the  ice  cover.  Thus,  a  fc*-type  value  that  corres¬ 
ponds  to  a  thin  floating  plate  will  be  much  larger  than 
the  k*  value  for  a  thick  plate  that  is  made  of  similar 
ice. 

The  post-buckling  branch  depends  on  the  magni¬ 
tude  of  the  fc*-type  value.  Therefore,  some  buckling 
loads  obtained  from  tests  with  very  thin  ice  plates 
may  not  be  suitable  for  the  interpretation  or  predic¬ 
tion  of  the  buckling  response  of  thick  ice  covers.  The 
results  shown  in  Figure  5  also  suggest  that,  when  using 
thin  plates  in  tests,  special  attention  should  be  devoted 
to  the  possible  occurrence  of  lift-off  and  its  effect  on 
the  obtained  results. 
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